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Abstract. Laplace's method is one of the fundamental techniques in the as- 
ymptotic approximation of integrals. The coefficients appearing in the re- 
sulting asymptotic expansion, arise as the coefficients of a convergent or as- 
ymptotic series of a function defined in an implicit form. Due to the tedious 
computation of these coefficients, most standard textbooks on asymptotic ap- 
proximations of integrals do not give explicit formulas for them. Nevertheless, 
we can find some more or less explicit representations for the coefHcients in the 
literature: Perron's formula gives them in terms of derivatives of an explicit 
function; Campbell, Froman and Walles simplified Perron's method by com- 
puting these derivatives using an explicit recurrence relation. The most recent 
contribution is due to Wojdylo, who rediscovered the Campbell, Froman and 
Walles formula and rewrote it in terms of partial ordinary Bell polynomials. 
The main purpose of this paper is to provide a new representation for the co- 
efficients, which contains ordinary potential polynomials. The proof is based 
on Perron's formula and Lagrange interpolation. The resulting formula has 
the advantage that these ordinary potential polynomials have a simpler gener- 
ating function than the partial ordinary Bell polynomials, which are required 
when using Wojdylo's formula. Due to the fact that the generating function 
is simple, it is possible, in many specific cases, to obtain an exact explicit for- 
mula for the ordinary potential polynomials, and therefore for the coefficients. 
The asymptotic expansions of the gamma function and the incomplete gamma 
function are given as illustrations. 



1. Introduction 

Laplace's method is one of the best-known techniques developing asymptotic 
approximation for integrals. The origins of the method date back to Pierre- 
Simon DE Laplace (1749 - 1827), who studied the estimation of integrals, arising 
in probability theory, of the form 

(1.1) IW^ e-^^^'''^gix)dx (A-^+oo). 



Here (a, b) is a real (finite or infinite) interval, A is a large positive parameter 
and the functions / and g are continuous. Laplace made the observation that the 
major contribution to the integral / (A) should come from the neighborhood of the 
point where / attains its smallest value. Observe that by subdividing the range of 
integration at the minima and maxima of /, and by reversing the sign of x whenever 
necessary, we may assume, without loss of generality, that / has only one minimum 
in [a, b] occuring at 2; = a. With certain assumptions on /, Laplace's result is 



f e-^-^(")<?(a;)da;-5(a)e-^/(")y^ 



2A/" (a) ■ 
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The sign ~ is used to mean that the quotient of the left-hand side by the right- 
hand side approaches 1 as A — > +00. This formula is now known as Laplace's 
approximation. A heuristic proof of this formula may proceed as follows. First, 
we replace / and g by the leading terms in their Taylor series expansions around 
X = a, and then we extend the integration limit to -foo. Hence, 

-'f^^^gix)dx^ re-K^('')+^(--'')^)5(a)da: 

J a 

n-\-00 ,1 

g[a)e ^f^^^ I e-^^^^~'^^' dx 
5(a)e-^^('^) 



2A/" [a] ■ 

The modern version of the method was formulated in 1956 by Arthur Erdelyi 
(1908 - 1977), who applied Watson's lemma to obtain a complete asymptotic ex- 
pansion for the integral (|l.ip . His method requires some assumptions on / and 
g. Suppose, again, that / has only one minimum in [a, b] which occurs at a: = a. 



(1.2) f{x)^f{a) + Y,akix~af+", 

k=0 

and 



(1.3) g{x) ^^bk{x - a)'' 



Jk V-^ ~ " 

with a > 0, 5R (/?) > 0, and that the expansion of / can be term wise differentiated, 
that is, 

00 



(1.4) /'(x) ~^afc(A:-f a)(a;-a)' 

k=0 

as a; — > a+. We may also assume, without loss of generality, that ao 7^ and 60 7^ 0. 
The following theorem is Erdelyi's formulation of Laplace's classical method. 

Theorem 1.1. For the integral 

I{\)= / e-^^(^)(?(a;)da;, 

J a 

we assume that 

(i) f{x) > f{a) for all x G (a, b), and for every d > the infimum of f{x) — f{a) 

in [a + (5, b) is positive; 
(a) f'{x) and g{x) are continuous in a neighborhood of x — a, expect possibly at 
a; 

(Hi) the expansions (|1.2p . (jl.3p and (jl.4p hold; and 
(iv) the integral / (A) converges absolutely for all sufficiently large A. 
Then 

('■^) nA)^e-""f;r(!l±i?)^. 

n— ^ ^ 

as A — > +00, where the coefficients c„ are expressible in terms of an and bn- 
The first three coefficients c„ are given explicitly by 

_ bo _ I (bi (/3 + l)ai&o 



(/3+i)/a a^ao 
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and 

To understand the origin of these coefRcients, we sketch the proof of Erdelyi's theo- 
rem. Detaild proofs are given in many standard textbooks on asymptotic analysis, 
e.g., in Erdelyi's original monograph [5, p. 38] or the classical books of Olver (TUl 
p. 81] and Wong [501 P- 58]. By conditions (ii) and (iii) there exists a number 
c e (a, h) such that f'{x) and g(x) are continuous in (a, c], and f'{x) is also positive 
in this interval. Let T = f{c) — f{a), and define the new variable t — t{x) as 

t^f{x)~f{a). 

Since f{x) is increasing in (a, c), we can write 

(1.6) /"'e-^/(")5(x)d2; = e-^-^(°) C e-^%[t)At 

J a Jo 

with h{t) being the continuous function in (0,T] given by 

(1.7) h{t)=g{x)—- 



fix)- 



By assumption, 



afe — a) as a; — >■ a . 

fc=0 



And so, by series reversion, we obtain 

oo 

x-ar^Y^ dkt''/" as t ^ 0+ 

k=l 

Substituting this into (|1.7p and using the asymptotic expansions (jl.3l) - (ll.4p yields 



k=0 



as t — >■ 0+. (In the common case where g = I, we have /3 — 1 and = 
(fc + 1) dfc+i/a.) We now apply Watson's lemma to the integral on the right-hand 
side of (|1.6p to obtain 

f e-^(^). (.) d. ^ e-/(^) g r (Ii±^) 

as A — -|-oo. To complete the proof, one shows that the integral on the remaining 
range (c, b) is negligible. 

The coefficients c„ are traditionally calculated manually in the way sketched 
above: by computing the reversion coefficients d„, substituting the resulting series 
into the asymptotic expansion of (jl.7L expanding the series in powers of t, and 
equating the coefficients in each specific application of Laplace's method. Due to 
this somewhat tedious computation, most standard textbooks on asymptotic ap- 
proximations of integrals do not give explicit formulas for the c„'s. Nevertheless, 
there are certain formulas in different degrees of explicitness for the coefficients 
Cn in the literature. When {x — a) " (/ (a;) — f (a)) and {x — a)^ ^ g {x) are an- 
alytic at a; = a, Perron's formula gives the coefficients in terms of derivatives of 
an explicit function involving / and g. Campbell, Froman and Walles rediscovered 
Perron's method and went further by computing these derivatives using an explicit 
recurrence formula The most recent contribution is given by Wojdylo, who 
rediscovered the Campbell, Froman and Walles formula and rewrote it in terms of 
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partial ordinary Bell polynomials |18| |19) . Using new ideas of combinatorial analy- 
sis, he was able to simplify and systematize the computation of the c„'s. We discuss 
these formulas in details in Section [2] 

Once we become familiar with the partial ordinary Bell polynomials, the deriva- 
tion of Wojdylo's formula from Perron's formula seems somewhat straightforward. 
In a bit trickier way, we obtain a new representation for the c„'s from Perron's 
formula which contains ordinary potential polynomials, in Section |31 The result- 
ing formula is similar to that of Wojdylo, but has an important advantage. The 
ordinary potential polynomials appearing in our new formula have a simpler gener- 
ating function than the partial ordinary Bell polynomials which are required when 
using Wojdylo's formula. Due to the fact that the generating function is simple, 
it is possible, in many specific cases, to obtain an explicit formula for the ordinary 
potential polynomials in terms of some elementary expressions or known combina- 
torial sequences. We will also give a representation for the partial ordinary Bell 
polynomials in terms of the ordinary potential polynomials. This representation 
shows that whenever one has an explicit expression for the latter type of polyno- 
mials, one has an explicit expression for the former type, even though it is more 
elaborate than the one for the latter type. 

In Section |4l we give two illustrative examples to demonstrate the application 
of our new method and to compare the results with those given by Wojdylo's 
formula. In the first example, we obtain several explicit expressions for the Stirling 
coefficients appearing in the asymptotic expansion of the gamma function. In the 
second example, we investigate certain polynomials related to the coefficients in the 
uniform asymptotic expansion of the incomplete gamma function. 

The definition and basic properties of the partial ordinary Bell polynomials and 
the ordinary potential polynomials are collected in Appendix El 



2. Perron's and Wojdylo's formula 

The asymptotic theory of integrals of type (|l.ip is also well established when A is 
complex and /, g are holomorphic functions in a domain of the complex plane con- 
taining the path of integration joining a to 6. A well-known method for obtaining 
asymptotic expansions for such integrals is the method of steepest descents. This 
method requires the deformation of into a specific path that passes through one 
or more saddle points of / such that the function 3 (/) is constant on it. (Recall 
that zo is a saddle point of / iff f {zq) — 0.) This new path is called the path 
of steepest descent. However, in many specific cases, the construction of such a 
path can be extremely complicated. This problem may be bypassed using Perron's 
method which - by requiring some extra assumptions ~ avoids the computation of 
the path of steepest descent, and provides an explicit expression for the coefficients 
in the resulting asymptotic series [Il][ini P- 103]. A direct adaptation of Erdelyi's 
theorem to complex integrals was formulated by Olver, where the coefficients in the 
asymptotic expansion can be computed formally in the same way as in the real case 
[TUl p. 121]. In addition, both Perron's and Olver's method allow the functions / 
and g to have algebraic singularities at the endpoint a with convergent expansions 
of the form dO]) and 

Based on Wojdylo's work, we show below that Perron's explicit formula also 
holds formally in the case of Erdelyi's theorem, when the expansions (|1.2p and 
(II. 3p may be merely asymptotic. Starting from the equation h (t) dt — g (x) dx, 
Wojdylo showed that 

{n+l3)/a n , .[,,,./ oo \ -(n+Z^)/"! 

dof % ' " a bn 

-C' 



bo fo ^" 




x=0 
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where d'^/dx'^ is the formal fcth derivative with respect to x. Through his analysis, 
he used these scaled coefficients c* . We shall not, however, use them in our paper. 
His expression can be simplified to 



(2.1) 



1 " 

{n+l3)/a X/ 



fc=0 



bn-k 
"fcT 



1 + E 



1 " 

{n+l3)/a X! 



fe=0 



k\ 



fe=l 



da;*' \Y,kLo^kx''+°' 



-{n+P)/a 



{n+l3)/a 



x=0 



aa. 



1 " 

n+B)/a E/ 



bn- 



("+/3)/ 




fe=0 



fc! 



J x=0 

{n+l3)/a 



If we identify / with its asymptotic expansion (|f .2^ , we arrive at the final forms 



(2.2) 



aa, 



1 " 

(n+,3)/Q X! 



6n- 







fc=0 



1 



,u(n+/^)/« 



d" 
dx" 



ao {x - af 



dx^ \f{x)-f{a) 

_ a \ (n+/3)/Q 

Gix) ' 



ao (x — aY 



fix) -Ho) 



where G {x) is the formal power series "YlV^obk {x — a)^ . An alternative form of 
this expression - using the formal residue operator - is stated in [iTJ p. 45]. If 
[x — a) " (/ (a;) — / (a)) and {x — a)^ ^ 9 {x) are analytic at x = a, this is the 
Perron formula. Using the definition of the ordinary potential polynomials in (|2.ip , 
we obtain 



1 " 

""o k=a 



(n+l3)/a,k 



ai 02 



1 5 ■ • • ; 



J /a; 



-Bfcj- (ai, 02, . . . , afe_j+i). 



This is essentially Wojdylo's formula. To make it more convenient to use, we apply 
the identity 



= (-1)' 



J + P - 1 
.1 



r{j + p) 



which yields 
(2.4) 



ar 



J_ v^, (-1)'^ Bfcj(ai,a2,...,afc_j+i) / n + /3 

Z^°"-'^'2^ («+/3)/a+i „■! la 

) k=a j=Q «o J \ 



+ j 



This is an even more exact form of Wojdylo's formula. In the common case, when 
5=1 (and therefore (3 — 1), formula (|2.4p is reduced to 



(-1) ^n.k (ai,a2, . . . ,a„_fc+i) f n+1 



(25) c ' V ^-'^ 



a ) k=0 "0 



fc! 



r 



Wojdylo's formula (|2.4I) . together with the recurrence (|A.ip for the partial ordinary 
Bell polynomials, provides a systematic way to calculate the coefficients c„. 
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3. New explicit formula for the c„'s 

Through this section, we identify / with its asymptotic expansion (jl.2p . Define 
the sequence of functions [Pk {z))j,^^ by 

Now, formula (j2.3D can be written as 

(3.1) „ J»+ff)/c ^^"-fe-P'^ 

Since by definition 



(n+f3)/c 
'■0 fc=0 



the Pk is a polynomial in z of degree at most k. When z = — j, j being a nonnegative 
integer, we can easily compute the value of these polynomials using the ordinary 
generating function 

( f{.x)-f{a) \ /^T , Ofc / NfcV 

7 Y^] =1 + 2^— 2^ -a ^2^Pk{-j){x-a). 

\a,{x-a) J y ^^ao ) t^. 

To obtain an explicit formula for the polynomial Pk in the whole complex plane, 
we use Lagrange interpolation. Taking the nodes to be {—j,Pk {—j)) (0 < j < k), 
we find 

Substituting z = (n + /3) /a and using the definition gives 

^ ^ (£1, - - 

Plugging this into p.ip yields the main result of this paper: 
Theorem 3.1. The coefficients Cn appearing in (|1.5p are given explicitly by 



^ ^k+3 n + 13 + ak fk 

w/iere Aj^fe = A^-^a; {a-ia^^ ,a2aQ^ , . . . ,afcaQ ^). 

In the common case when g = 1, we have P — 1 and formula p.2p is reduced to 
r{n±L+n+l) " (-if /n^ 



(n+l)/"„| r f!l± 



auQ n 



I Q / fc=o Q 



(3.3) 

(n+l)/a\ ^ ) A^^ „ _|_ 1 _|_ Q,fc \^ 
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From the recurrence relation (|A.2p for the ordinary potential polynomials it follows 
that 

n 
1=0 

This, together with p.2p . provides a simple and efficient way to compute the coef- 
ficients c„. 

In many special cases it is possible to obtain an explicit formula for the ordinary 
potential polynomials (or the partial ordinary Bell polynomials) in terms of some 
elementary expressions or known combinatorial sequences. 

Before we proceed with examples, we derive connection formulas between the 
quantities Bkj = B^j (01,02, . . .,ak-j+i) andAj^k = Aj,fc {aia^^ ,a2ao^ , . . . ,0^00 ^) 
appearing in Wojdylo's and our formula, respectively. The corresponding generat- 
ing functions are 

(00 
k=l 

and 

(00 

Since the Aj-.^'s have simpler generating function, it seems plausible to express the 
partial ordinary Bell polynomials ^ in terms of the ordinary potential polynomi- 
als Aj.fc. By definition 

A.. - E 01 B.. ^, . . . , ^) = ± {^) W. 

(This is true also when k < j since B^ ^ = for k < i < j.) This means that the 
sequence {Aj^k)j is the binomial transform of the sequence ^(—1)"' OQ-'Bfe . Since 
the binomial transformation is an involution, we have the required expression 

(3.6) Bu,, = {-iyaij2{-ir(^\K,k. 

Naturally, this can also be deduced from the generating functions p.4p and p.5p 
using the binomial identity. 

4. Examples 

In this section, we give two illustrative examples to demonstrate the application 
of Theorem 13.11 as well as to compare the results with that of given by Wojdylo's 
formula. In the first example, we derive explicit expressions for the so-called Stirling 
coefficients appering in the asymptotic expansion of the gamma function. In the 
second example, we obtain an explicit formula for certain polynomials related to the 
coefficients in the uniform asymptotic expansion of the incomplete gamma function. 

Example 1. The gamma function can be defined by the following integral 

p + OC 

(4.1) r{X + l)^ e-H^dt, A>0. 

^0 

If we put t — X(l + x), we obtain 

/+00 
g-A(:.-log(l+:r))^^^ 



J ^°" Uo(x-a)° 

.kX f fix)-fia) Y ^ , 
-a = 7 = Aj- fc X - a . 
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and hence, using the identity F {X + 1) ~ XF (A), 



(4.2) 



A^e--^ 



+ 00 



foUows. Consider the first integral Let / (x) — x — log(l + x), x >Q. Then / has 
a global minimum at a: = and the expansion 



/(x) = ^(-l) 



k X 



k+2 



fe=0 



k + 2 



holds as X — ^ 0+. We can apply Theorem 11.11 with g = 1, a = 2, /3 = 1, bo = 1, 
&fc = for fc > 0, and ak (-l)*" / {k + 2). The resuh is 



(4.3) 



p-\-oo / _i_ 1 \ 

-^0 n \ / 



A("+i)/2 



where, by Perron's formula. 



(4.4) c, 
Similarly, one finds that 



1 
2n^. 



dx" \x — log (1 + x) 



(n+l)/2' 



where 



2n\ 



d" 
da;*^ 



-X ~ log (1 — x) 



ri=0 



(n+l)/2' 



n+1 



A("+i)/2 



til 

2nl 



J x=0 
^2 



da;" \x ~ log (1 + a; 



By substituting these asymptotic series into 
arrangement, we deduce 



(n+l)/2' 

= (-!)" c„. 

- x=0 

and performing a simple re- 



(4.5) 

as A — > +0O. Here 
(4.6) 



r(A)^V2^A^-V2e-A^(_l)«^ 



"+2 



C2r, 



2"n! 



n=0 



]2n 



da;2" \2a;-log(l + a;) 



n+l/2' 



are the so-called Stirling coefficients [T31 p. 26]. This representation is well-known 
(see, e.g., [I] [23 P- HI])- The first few are 70 = 1 and 

1 1 139 571 



71 



72 



73 



74 



12' ' 288 ' 51840 ' 2488320 
In this case, we are able to obtain explicit formulas for both B„.fc and Afc.„. Due 
to the simpler generating function, we use our new formula (j3.3p and compute the 
ordinary potential polynomials. Here the generating function p.5p takes the form 

k 00 



X - log (1 + x) 
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where kk,n = kk,n (aiOQ \a2ao \ • ■ • , a„aQ ^) and a„ = (-1)" / {n + 2). Using the 
generating function of the (signless) Stirhng numbers of the first kind yields 

(x - log (1 + x))' - (^) (- (1 + ^))' 



3=0 
k 



V 7 / 

j=0 ^-^ ^ »i=0 



-k+j fk\.,s{n-k + j,j) 

n=0 \i=0 



j/ ' {n-k + jy. 



which gives 

Substituting this into (I3.3P produces 

2"-V2r(3n+|) ^ 2^- /2n\^ - fk\ ., s{2n + k + j,, 

2n 



_ ^ 2"+^+i/^r(3n+§) (-1)^- s (2n + fc + j, j) 



r (n + i) ^ (2n + 2fc + 1) (2n - fc)! ^ (fc - j)! (2n + fc + j)! ' 
Finally, by g^l), 

(-ir2"+^-+ir(3n+|) ^ („i)^-,(2n + fc + ,,j) 
^ ■ "^'^ ^ V^(2n + 2/c + l)(2n-/c)!^^ (fc-j)!(2n + fc + j)! ■ 

As far as we know, this representation of the Stirling coefficients is entirely new. 
Using expression (I3.6P together with (j4.7p . Wojdylo's formula (j2.5p yields the more 
complicated formula 

^^ (-ir2"r(n + fc+l) ^ 2^- ^ (-l)'.(2r^ + J+^,^) 
^ ■ ^ ^" to ^ ^^(fc-j)!^ (j-*)!(2n + j + *)!' 

which is the result of Lopez, Pagola and Perez Sinusfa 

We remark that the substitution x = log(f/A) in (|4.ip leads to the form 

/^°°e"^«=^--'i)dx+ r°°e-^(°"^+--i)dx. 
e Jo Jo 

Similar procedures to the one described above give (|4.5I) with 

(-ir 



7n 



2"n! 



d2» /I x2 ^ 



dx^" \ 2 e^ — X — 1 

J a;=0 



' t 0F(2n + 2fc + l)(2n-A:)! t (/c- j)!(2n + /c + j)! 

_ 1^ (-1)" 2-r (n + fc + 1) ^ 2^- ^ (-1)' ^ (2n + j + z) 
^ ■ ^ t^ ^t (fc-j)! t 0-~z)!(2n + j+z)! ' 

using formula (I2.2p . p. 31) and (12.51) . respectively. Here S{n,k) denotes the Stirling 
numbers of the second kind. The first representation is known (see, e.g., [T] [3]). 
The second formula is a simplified form of the one derived in [5]. To our knowledge. 
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the third one is entirely new. Note the remarkable similarity between expressions 
(gSl) and and expressions and (|ilT|) . 

Example 2. The incomplete gamma function is traditionally defined by the fol- 
lowing integral: 



r{a,x)= / e'H"-~'^dt, a > 0, x > 0. 

J X 

For large values of x, the function admits the following asymptotic expansion 

. T-n T -, a—1 (a— l)(a — 2) 
r (a, a;) "e'' - 1 + + ^ ^ + ■■■ , 

X X'^ 

which is useful only when a — o{x) 11, p. 179]. However, in exponentially improved 
asymptotics, we need the asymptotic properties of F (a, a;) as a — ?> +00 and x — Xa, 
where A ^ is a constant. (Do not confuse it with the variable of the integral 
(|l.ip .) Starting from the integral representation 

r (n r) p-'^v(A) /-c+ioo 1 

= / e-^-W-i-dt, < c < A, 

r (a) 27ri J^_i^ A - t 

^(t) -t-logi-1, 

Temme [16j gave the following uniform asymptotic expansion of the (normalized) 
incomplete gamma function as a — +00 

r{a,x) 1 / Ar\ e-^ ^ C„ (r;) 



r{a) 2 V V 2 / ^r, a" 

^ ^ \ / n— 



Here 77 = ^y2lp{X) and erfc denotes the complementary error function. He gave 
a recurrence for the coefficients C„ (77) and showed that they have the general 
structure 



(4.12) C„ (r;) = (-1)' 



)„(^) 2»r(n+i)~ 



where /i = A — 1 and (5„ is a polynomial in /i of degree 2n [T7] . The first few of the 
Qn's are given by 

Qo (m) = 1: 



Qi (m) = 1 + M+ 

25 1 1 

Q2(/.)=3 + 5m+^M^ + ^M^ + ^//. 



As far as we know, no simple explicit formula for these polynomials exists in the 
literature. Dunster et al. W derived the integral representation 

(4.13) a. <" + «^ 



(27r)^/^ i (z- A) (z - log z - 1)"+i/2 ' 

Here, is a loop in the z plane enclosing the poles z = 1 and z = X. The term 
{z — logz — i)"+i/2 ^^^Y and positive for argz = and z > 1 and is defined by 
continuity elsewhere. We suppose that A ^ 1 (77 7^ 0) and use the Residue Theorem 
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to obtain 



Cn iv) = 
= (-1)" 
= (-1)" 



/2tt 

r(n + i) 



Res 

2 = 1 



1 



1 



J 2)1 



{z - A) (z - log z - 1)"+^/' (A - log A - 1)"+^/' 

2 \ n+1/2 ' 



1 



(z-ir 



1 



22n+l/2^| 



dz^" j (z — A) \ z — log z — 1 
1 



2"r(»+^) 



dcc^"- I {^i — x) \a; — log (1 + a;) 



n+l/2~ 



z = l 



J an=0 



,2n+l 



2n+l 



Upon comparing this with (|4.13p . we see that 



2n+l 



22n+l/2„I 

and from Perron's formula 



J2n 



in - x) \x - log (1 + x) 



1+1/2^ 



x=0 



(4.14) 



(2n)!/i' 



|,,2n+l 



22n-l/2„! 



■C2n 



r n 



,,2ji+l„ 



where the c„'s are the coefficients appearing in the asymptotic expansion of the 
integral (jl.ip with a — 0, b > 0, f {x) = x — log (1 + x) and 5 (x) = {jj, — x)^^. 
Now, using Theorem 13.11 we derive an explicit formula for the polynomials Qn in 
terms of the Stirling numbers of the first kind. Near x — g {x) — X]fc°=o 
therefore bk — /i^*^^^. The corresponding ordinary potential polynomials are given 
by (|4J| . Hence, by Theorem [Q 



C2r, 



1 ^ 2"+l/2r (f 

2/ fc=0 

This, together with (|4.14l) . yields 



rin 



kl 



■E 



i-iy 



2n + 2j + 1 \ j 



2n 



2"+ir I 



fc=0 
^2n An) k 



0Ffc! 



E 



(-1)^ 



J=0 



2n + 2j + 1 Vj 



or Qn ifJ.) = Efe=o ^fe V", where, by dUT]), 
(4.15) 



(n) ^ (-1) 2"+^+ir (n + fc + I) ^ (_!)' g (fc + j + i) 



/7r(2n + 2j + 1) (A:- j) 



(j-i)!(/c + j + i)! 



We remark that = (—1)" 7n, where 7„ is the nth Stirling coefficient. Naturally, 
we could use Wojdylo's formula (|2.4p to obtain an aternative representation for the 
coefficients c„, and thus for the polynomials Qn (/i) themselves. The corresponding 
partial ordinary Bell polynomials can be given explicitly using expression (|3.6p 
together with (14. 7p . the resulting formula is, however, even more elaborate than 
(SIB. 

When A = 1 (77 = 0), we have C„ (0) = lim,;^o C-'n (??) and Tcmme's expansion 
takes the form 

r{a,a) _ 1 1 ^ Cn (0) 
2 2^ 



r(a) 



/27ra ^ a" 

n—O 



as a — > +CX}. However, the value of the limit lim,j_>.o Cn (v) is not obvious from 
(|4.12p . To derive an explicit formula for the C„ (0)'s, we proceed as follows: for 
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every positive integer m 

r (m, m) „, ^^^—^ TO™e~™ TO™e~™ / m\ k\ 

— !^ — . — '- = e~ > = \ > ' ' 

r (m) ^ kl ml ml ^ 

^ ' fc=0 fc=0 



m\ ml ^-^ \ k 

k=0 



+00 



+00 



It is well-known that 



(4.16) '-^^ T^E— ' 

ml ^/2^^m" 

^ n— 

where 7„ denotes the nth Stirhng coefficient [131 P- 26]. By (14. 3p . we have 



-l + m e-"(^"'°s(i+^»da; - -1 + V T 
•^0 n=0 



7"/2 ' 



where c„ is given by (|4.4p . From these we deduce 

r-'/ \ -.00 00 oo/ _, 

r (m, to) JiL _| 3_ y^7ri_y^^/'^+l 



^/n-2fc+l\ \ 1 



^ TO" ^ 1 ' V 2 y I TO"/2 



00 00 / Ln/2_ 

2.TO^„TO" ^^^y^^ 

1 1 °° 7 00 /L(«+1)/2J 

2 " V2^ ^ ^ V2^ \ S ^"'^ ^ 



fc + 1 ) C„_2fe+1 



l"/2 • 



The well-known identity J2'k=o ^ Ikln-k = (n > 1) [TH p. 33] gives 



I,— n ^ ' t — n 



fc=0 ^ ^ ' k=0 

from which we obtain 



r(m,m) 1 1 ^ / ^ ^, 

r(TO) 2 V2^;^oV .t^o 



n - /C H- 1) C2„-2fc+l ) 



TO" 



From the uniqueness theorem on asymptotic series we finally have 

n 

Cn (0) -7„ + ^ 7^/" (n - fc + 1) C2„-2fe+l 



n-1 



oTn + E ~ fc)!7fcC2„-2fe-|-l- 



k=Q 

We remark that one can start with the representation 

r(TO,TO) _ TO"'e-"' r°°g-™(.-log(l+.))^d^^ 

r (to) r (to) 7o 1 -I- x ' 

which follows from the definition of the incomplete gamma function. The factor 
before the integral can be expanded into an asymptotic series using (|4.16p . The 
asymptotic expansion of the integral can be deduced from Erdelyi's theorem. In 
this way, we get a similar representation for the coefficients C„ (0) to the one above. 
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However, in this case g (x) = {1 + x) , whereas in the previous one g = 1, which 
leads to a somewhat comphcated representation. 

5. Conclusion 

Laplace's method is one of the classical techniques in the theory of asymptotic 
expansion of real integrals. The coefficients c„ appearing in the resulting asymp- 
totic expansion, arise as the coefficients of a convergent or asymptotic series of a 
function defined in an implicit form. Traditionally, series inversion and composition 
have been used to compute these coefficients, which can be extremely complicated. 
Nevertheless, there are certain formulas of varying degrees of explicitness for c„'s 
in the literature. One of them is Perron's formula which gives the coefficients in 
terms of derivatives of an explicit function. The most explicit formula is given by 
Wojdylo, who rewrote Perron's formula in terms of some combinatorial objects, 
called the partial ordinary Bell polynomials. 

In this paper we have given an alternative way for simplifying Perron's formula. 
The new representation involves the ordinary potential polynomials which have 
simpler generating function than the partial ordinary Bell polynomials appearing 
in Wojdylo's method. Due to this simple generating function, in many specific 
cases it is possible to obtain an explicit formula for the coefficients c„ in terms 
of some elementary expressions or known combinatorial objects. Alternatively, the 
recurrence (|3.2p can be used for computing the ordinary potential polynomials, and 
hence the coefficients c„ themselves. 

We have applied the method to the two important examples of the gamma func- 
tion and the incomplete gamma function. We have obtained new and explicit 
formulas for the so-called Stirling coefficients appearing in the asymptotic expan- 
sion of the gamma function. We have also derived explicit formulas for certain 
polynomials related to the coefficients in the uniform asymptotic expansion of the 
incomplete gamma function. 

It turns out that formally the same coefficients appear in the asymptotic expan- 
sion of contour integrals with a complex parameter. Our new method, hence, can 
be applied to give explicit formulas for the coefficients arising in the saddle point 
method [TOl p. 125] [11, p. 47], the method of steepest descents [7] and even in the 
principle of stationary phase [5] [HI p. 45]. 

Appendix A. Powers of Power Series 

Let F {x) = 1 + fnx" be a formal power series. For any nonnegative 

integer fc, we define the partial ordinary Bell polynomials Bn^k (/i, /2, • • • , fn-k+i) 
(associated to F) by the generating function 

Coo \ oo 

E = E B".'^ (/l' • • ■ , /n-fe+l) a;", 

ri=l / n=k 

so that Bo,o = 1, B„,o = (n > 1), B„_i = /„ and B„,2 = YTjZl fjfn-j- From the 
simple identity, [F (x) - 1)''+^ = (F [x) - 1) (F {x) - if ' , we obtain the recurrence 
relation 

n—k 

(A.l) B„^fe+i (/l, /2, . . . , fn-k) = ^ fjBn-j.k (/l, /2, • ■ • , /n-j-fc+l)- 

J = l 

An explicit representation, that follows from the definition, is given by 

u 

D ( f f f \ — \ ^ f'^l f ^2 J-Kn-fc + l 

^n,k {J1,J2,---, Jn-k+1) — 7 , \i\ J2 ' ' ' Jn-k+1 ' 
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where the sum runs over all sequences /ci , ^2 , . . . , A;„ of non-negative integers such 
that ki + 2k2 + ■ ■ ■ + {n — k + I) kn-k+i = n and fci + fc2 + • • • + kn-k+i ~ k. 

For any complex number p, we define the ordinary potential polynomials Ap^„ (/i , /2, 
(associated to F) by the generating function 

(oo \ ^ oo 

n=l / ri=0 

hence, 



Ap^n (/l, /2, ■ • • , /n) — ^ L j B„^fc (/l, /2, ■ • • , /ri-fc+l)- 

The first few are Ap_o = 1, Ap_i = pfi, Ap_2 = p/2 + (2)/! ' ^^'i general 
Ap,„ (/i, /2, . . . , /„) - ^ k) fc^i^^TTTfcj/i'Vs^ ■ ■ • /n 



/ n ' 



where the sum extending over all sequences fci , fc2 , . . . , fc„ of non- negative integers 
such that ki + 2k2 + ■ ■ • + nfc„ = n and ki + k2 + ■ ■ ■ + kn — k. Since {F {x)Y = 
{F {x)Y ^ F (x), we have the recurrence 
(A.2) 

n 

Ap,n (/l, /2, ■ ■ ■ , fn) — ^p-l.n (/l, /2, • ■ • , /n) + ^ fkAp-l,n-k (/l, /2, • • • , fn-k)- 

k=l 

In general, if G (a;) = J2'^=o 5"^" is a formal power series, then by definition, we 
have 



00 / n 



G {y {F [x) - 1)) = E E 5feB„,fe (/i, /2, • ■ • , /„-fc+i) 2/'= x". 



Specially, 



n=0 \/c=0 



00 / n 



exp ( y E ) = E ( E ^-ii/llA^pA^.'c ,^ 

V n=l / n=0 \fe=0 ■ / 

For more details see, e.g., Riordan's book [TSl p. 189]. 
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